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FOREWORD 


The  objective  of  this  project  is  the  understanding  of  the 
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1.  INTRODUCTION 


In  recent  years  high-energy  propellants  have  been  developed 
containing  cyclotetramethylenetranitramine  (HMX)  to  improve  the 
performance  of  rocket  motors.  In  testing  such  propellants, 
several  accidental  explosions  have  occurred  during  rocket  motor 
firings .  The  goal  of  this  study  is  the  identification  of  the 
causes  of  the  explosions  in  order  to  prevent  such  occurrences 
in  the  future. 

This  program  examines  possibilities  of  deflagration  to 
detonation  transitions  (DDT)  caused  by  sudden  exposure  of  cracks 
to  high  pressure/ temperature  gas  cavities  such  as  the  combustion 
chamber.  It  consists  of  the  following  steps: 

•  development  of  analytical  means  with  which  to 
identify  crack  conditions  and  propellant  prop¬ 
erties  needed  to  produce  pronounced/rapid 
pressure  rises  within  burning  cracks 

•  comparison  of  predicted  pressure  transients 
with  stepwise  shock  wave  amplitudes  and 
durations  known  to  initiate  propellants ,  and 

•  validation  of  analytical  predictions  by  ex¬ 
perimental  means . 

Results  of  previous  IIT  Research  Institute  (IITRI)  studies  are 
described  in  references  [1]  through  [4]. 

1.1  Background 

During  previous  IITRI  studies  two  factors  were  identified 
as  possible  causes  of  pronounced  pressure  rises.  Namely  stress 
waves  that  partially  collapse  burning  cracks,  and  foam  or  melt 
layers  that  are  rapidly  consumed  by  dynamic  burning.  Stress 
waves  are  generated  by  rapid  pressure  rises  within  cracks.  They 
may  be  produced  by  sudden  exposure  of  cracks  to  high-pressure 
gas  cavities  and/or  rapid  burning  of  the  propellant  within  the 
cracks.  The  consequence  of  applying  stress  waves  to  a  burning 
crack  is  stress  waves  of  appreciably  greater  amplitude.  Resul¬ 
tant  stress  waves  will  then  compress  any  neighboring  burning 
cracks  and  produce  stress  waves  of  still  higher  amplitude.  By 


this  process,  progressively  more  pronounced  pressures  can  be 
generated  within  neighboring  burning  cracks  on  a  crack  by  crack 
basis . 

Pressure  rises  are  also  highly  dependent  upon  the  foam  mass 
and  temperature,  and  upon  the  size  and  orientation  of  the  cracks 
with  respect  to  each  other.  With  foam  masses  of  the  order  of 
0.01  g/cm^  or  greater,  only  a  few  parallel  closely  spaced  burning 
cracks  are  needed  to  achieve  pressures  of  the  order  of  10  kbars. 
Much  of  the  pressure  rise  occurs  within  a  fraction  of  a  ysec. 

1.2  Status  of  Present  Study 

The  scope  of  work  for  the  present  study  consisted  of  the 
following  analytical  and  experimental  tasks: 

(1)  Measurement  of  HMX  foam  mass  present  during 
steady  burning  at  ambient  pressure 

(2)  DDT  experiments  in  which  burning  propellants 
are  impacted  by  a  flyer  plate 

(3)  Incorporation  of  pressure-dependent  propellant 
impedance  in  computerized  models 

(4)  Prediction  of  transient  pressures  produced 
sequentially  within  a  series  of  burning  cracks 

(5)  Development  of  a  computer  code  with  which  to 
predict  amounts  of  foam  or  melt  generated  within 
cracks  suddenly  exposed  to  cavities  containing 
high  temperature/pressure  gases 

No  sources  of  HMX  cylinders  with  dimensions  of  the  order  of 
several  centimeters  could  be  found.  Thus  20  cylinders  of  pro¬ 
pellant  composed  of  approximately  65  percent  HMX  and  35  percent 
polyethlylene  glycol  binder  were  ordered  in  May  1979  for  the 
foam  mass  and  DDT  experiments  cited  by  items  1  and  2.  Unfor¬ 
tunately  the  propellant  cylinders  were  not  received  in  time  to 
conduct  the  experiments.  Instead  efforts  were  expended  in: 

•  upgrading  the  experimental  design  described  in  the 
previous  interim  report  [4] , 

•  assessing  the  effect  of  sudden  pressure  relief 
upon  foam  masses  present  during  steady  burning, 

•  assessing  rates  of  temperature  rise  of  gases 
evolved  by  heated  HMX  prior  to  ignition,  and 

•  analyzing  means  for  predicting  the  ignition  of  HMX. 
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Items  3  through  5  have  been  completed  and  resu?.ted  in  two  com¬ 
puterized  models  called  PROS  and  PROP.  PROS  is  a  new  code  that 
predicts  dynamic  burning,  melt  formation,  crack  deformation  and 
gas  pressures  along  the  length  of  cracks  following  sudden  ex¬ 
posure  to  high-pressure  gas  cavities.  Figure  1(a)  illustrates  the 
problem  treated  by  PROS.  Gases  initially  within  the  crack  are 
at  a  state  of  rest  with  a  uniform  pressure  and  temperature. 

Initial  crack  widths  may  be  uniform  or  variable  along  the  length 
of  the  crack.  Gas  flows  within  the  crack  are  treated  as  one¬ 
dimensional.  Gas  parameters  and  burning  conditions,  will  of 
course,  vary  with  crack  location  and  time  after  the  crack  propa¬ 
gates  into  the  cavity. 

PROP  is  a  revised  edition  of  the  code  used  in  the  previous 
interim  report  [4].  It  predicts  dynamic  burning,  crack  deforma¬ 
tions  and  gas  pressures  produced  within  an  element  of  a  burning 
crack  subjected  to  stress  waves.  Amounts  of  gas  escaping  from 
the  crack  element  are  considered  insignificant  over  the  times 
of  interest.  During  the  past  year  a  subroutine  was  incorporated 
into  this  model  in  order  to  guide  development  of  the  DDT  fixture. 
The  subroutine  includes  provisions  for  accounting  for: 

•  generation  of  a  system  of  plane  stress  waves 
radiating  from  the  burning  crack, 

•  reflection  of  the  radiating  stress  waves  by  a 
high-density  media  such  as  the  motor  case, 

•  partial  compression  of  the  crack  by  the  re¬ 
turning  stress  waves,  and 

•  subsequent  expansion  of  the  crack  by  ele¬ 
vated  gas  pressures  produced  within  the  crack. 

The  remainder  of  this  report  describes  basic  analytical 
techniques  used  by  the  two  mod-  Is,  the  revised  DDT  fixture,  and 
various  predictions  derived  from  the  analytical  studies.  Nomen¬ 
clature  and  property  data  are  presented  in  the  appendix  at  the 
rear  of  this  report. 
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Figure  1(a).  Crack  Configuration  Treated  by  PROS 


Figure  1(b).  Crack  Configuration  Treated  by  PROP 
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(Open  or  Closed) 


2.  BASIC  EQUATIONS  AND  ANALYTICAL  RESULTS 

This  section  briefly  describes  means  by  which  various 
phenomena  are  treated.  Also  included  are  a: 

•  discussion  of  ignition  criteria,  and 

•  predictions  of  changes  in  the  foam  mass  of 
steadily  burning  HMX  following  sudden  pres¬ 
sure  relief. 

2.1  Propellant  Heating 

Rates  of  propellant  burning  depend  upon  the  temperature  pro 
file  within  the  propellant,  which  in  turn  depends  upon  the  time- 
dependent  heat  fluxes  supplied  the  propellant  before  and  following 
the  start  of  burning.  Systems  of  equations  [4]  are  used  to  pre¬ 
dict  dynamic  burning  in  terms  of  the  heat  fluxes  received  by  the 
propellant  crack  surfaces. 

This  section  is  concerned  with  describing  heat-flux  expres¬ 
sions  developed  for  PROS.  In  this  regard,  the  equations  used  by 
PROP  are  special  cases  of  the  equations  used  by  PROS.  Heat 
fluxes,  will  of  course,  vary  with  the  physical  condition  cf  the 
propellant  surface,  i.e. ,  solid,  melting  or  burning. 

2.1.1  Convective  Heat  Flux  q  Prior  to  Start  of  Melting 

Prior  to  the  start  of  melting,  the  incident  convective 
heat-flux  q  is  described  by  the  expression  used  by  ICuo  [5].  In 
terms  of  the  gas  temperature  T^  and  the  temperature  T^  of  the 
propellants  surface,  the  heat  flux  q  is 

q  -  VVTf> 

where  the  heat  transfer  coefficient  hc  is  given  by 

n  i  ns 


(1) 


The  dynamic  heat  flux  q  described  by  equations  (1)  and  (2)  is 

used  in  PRO'’.  A  constant  heat  flux  cs  is  used  in  PROP. 

’  ..  .  '  ■) 

u  2.1.2  Heat  Flux  q^  During  Melting 

Two  sources  of  heat  flux  are  present-  while  the  propellant 
surface  is  melting.  The  first  is  the  heat  flux  q  described  in 
subsection  2.1.1;  the  second  is  heat  generation  within  the  melt 
or  foam  layer.  The  latter  is  termed  internal  heating.  While 
the  propellant  is  melting  prior  to  burning,  the  heat  flux  q^  is 
described  by 


U  ' 


u 


c 
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q£  =  q  +  P  rf  Qs 


(3) 


The  term  q  represents  the  incident  heat  flux,  while  the  second 
expression  represents  internal  heating.  In  that  the  incident 
heat  flux  is  convective  it  will  decrease  as  the  foam  temperature 
rises.  On  the  other  hand,  increased  foam  temperature  accen¬ 
tuates  internal  heating  by  raising  the  propellants  regression 
late  r^.  The  latter  points  up  the  importance  of  internal  heating 
in  speeding  the  melting  of  the  propellant  and  in  accelerating 
subsequent  burning. 


2.1.3  Heat  Flux  q^  During  Burning 

Gases  evolved  by  burning  propellants  interact  with  gtses 
flowing  laterally  along  crack  surfaces,  and  thereby  perturb  the 
convective  heat  fluxes  cited  earlier.  Heat  fluxes  resulting 
from  the  interaction  of  the  two  gas  flows  are  termed  erosive. 

In  this  regard,  Lenoir  and  Robillard  [6]  developed  an  expression 
for  the  rate  r^  of  steady  burning  produced  by  erosive  heating. 
Subsequently  their  expression  was  used  by  Kuo  [5].  It  is  given 
by 


n  Bprf 

rf  =  a  P  +  k  hc  exp(-~-~  ) 

E 


(4) 


Since  then  King  [7]  upgraded  the  first  express  i  on  the  right- 
hand  side  of  equation  (4)  to  account  for  the  effect  of  erosive 
burning  upon  the  flame  standoff  distance.  The  result  is 
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faPn12  Bprf 

rf  *  ^-  +  -*hcWt-  -inf)' 

rf  s 


During  steady  burning,  the  heat  flux  must  maintain  the  . 
foam  temperature  at  a  steady  value  so  that 

5f=  (6) 


For  purposes  of  brevity,  the  expression  within  the  brackets  of 
equation  (6)  shall  be  represented  by  Qf.  Multiplying  equation  (5) 
by  p($£  and  replacing  P?fQf  by  its  equivalent  q^  yields 


n\ 2tt 


P  (aPn)  Q.c 


- - "  +  khcexp  (-  -rF)P^f 


When  the  cross-flow  gas  velocity  u  is  zero,  the  second  expression 
on  the  right-hand  side  of  equation  (7)  vanishes  and  the  steady 
heat  flux  qf  equals  the  first  expression  given  by  equation  (7). 

In  order  to  describe  the  heat  flux  qf  present  during  dynamic 
burning,  we  shall  replace  the  first  expression  on  the  right-hand 
side  of  equation  (7)  by  the  expression  developed  by  Krier  [8]  for 
dynamic  burning  as  follows: 


P(aPn)2Qf 


p (aPn)2 (Qf  “  Qg) 


+  prf(prfQg) 


pr-(C  -  C  )  (T-  -  T-) 
y  f  m  g  f  f 


Notice  that  the  above  expressions  are  identical  for  steady 
burning  wherein  r^  =  r^  =  a?n  and  Tf  =  T^. 

According  to  Lenoir  and  Robillard  [6]  the  factor  k  given 
in  equation  (7)  equals 


Substituting  the  expressions  given  by  equations  (8)  and  (9)  into 
equation  (7) ,  and  replacing  q^  and  r^  by  their  unsteady  values 
q£  and  r^  yields 


p(aPn)2«}f-Q  ) 

qf - 1  '  +  PrfQs  +  Prf(Cm-  Cg)(Tf  -Tf)  + 


6prf 

hc(WeXP(-  M>  (10) 

The  dimensionless  constant  8  of  equation  (10)  must  be  evaluated 
experimentally.  In  that  its  value  for  HMX  is  not  known  to  the 
author  it  shall  tentatively  be  assumed  equal  to  the  value  15 
determined  by  King  [7]  for  a  propellant  different  from  HMX. 

2 . 2  Dynamic  Flow  of  Gases  Within  Cracks 

An  essential  aspect  of  the  problem  is  the  prediction  of  the 
transfer  of  gases  and  heat  along  the  length  of  a  crack.  For 
this  purpose  a  gas  dynamic  model  was  developed  [1].  It  is  based 
upon  an  Eulerian  representation  of  the  flow  region  and  is  limited 
to  a  one- dimensional  treatment.  The  Eulerian  representation  was 
selected  in  that  it  permits  simple  identification  of  crack  loca¬ 
tion  and  allows  for  mixing  of  freshly  burnt  reaction  gas. 

One-dimensional  treatment  is  appropriate  in  that  crack 
thicknesses  and  lengths  are  usually  large  compared  to  their 
width.  Thus  the  most  significant  flow  gradients  and  movements 
occur  along  the  length  of  the  crack.  Values  used  in  the  one¬ 
dimensional  treatment  represent  cross-sectional  values. 

The  model  treats  a  crack  with  a  specified  width  variation 
at  zero  time.  The  gas  in  the  crack  at  zero  time  is  considered 
to  have  the  same  composition  as  the  reaction  products,  and  to 


exist  in  a  rest  state  at  some  uniform  temperature  and  pressure. 
The  crack  is  then  suddenly  connected  to  a  cavity  filled  with 
gaseous  reaction  products  at  high  temperature  and  pressure.  The 
cavity  state  is  treated  as  constant  with  respect  to  time  due  to 
the  relatively  short  duration  of  the  events  occurring  within  the 
crack. 

The  downstream  end  of  the  crack  can  be  either  closed  or 
connected  to  another  cavity  held  at  a  specified  pressure.  In 
either  case  crack  length  is  held  constant. 

2 . 3  Mechanical  Deformation  of  Cracks 

There  are  two  sources  of  mechanical  response  of  cracks.  The 
first  is  elevated  gas  pressures  within  cracks  that  act  to  expand 
cracks.  The  second  is  due  to  stress  waves  radiating  from  cracks 
and  subsequently  interacting  with  other  cracks  or  the  same  crack 
following  interaction  of  the  radiating  stress  waves  with  the 
confining  shell  structure. 

The  mechanical  response  of  the  propellant  mass  to  pressures 
within  cracks  is  very  complex  due  to  the  multiplicity  of  stress 
wave  types  and  the  complex  geometry  of  the  crack,  propellant 
mass  and  confining  shell.  The  primary  motion  of  the  propellant 
mass  will  be  in  a  direction  normal  to  the  crack  and  will  be 
roughly  one  of  plane  strain.  For  this  reason  a  simple  local 
plane  strain  representation  for  the  transient  stress  field  was 
adopted  in  the  absence  of  remote  boundary  effects. 

A  simple  shell  response  model  was  selected  to  treat  the 
remote  boundary  effect  upon  reflected  stress  waves.  The  model 
treats  the  shell  as  a  simple  single  degree-of- freedom  lumped 
mass  which  is  position  restrained  and  which  responds  to  the 
mean  radiative  stress  field.  The  response  generates  a  system 
of  reflected  stress  waves  which  arrive  back  at  the  crack  after 
a  delay  period  which  depends  upon  the  distances  to  the  remote 
boundary . 


Equations  describing  the  velocities  of  the  two  crack  walls 
caused  by  gas  pressures  and  incident  stress  waves  are  presented 
in  the  previous  interim  report  [4] . 

2.4  Ignition  Criteria  and  Foam  Mass 

2.4.1  Ignition  Criteria 

Ignition  refers  to  the  start  of  propellant  burning.  Knowl¬ 
edge  of  when  ignition  occurs  is  important  in  that  the  heat  flux 
undergoes  rapid  change  following  ignition. 

Ignition  represents  a  complex  phenomenon  which  we  attribute 
to  exponential  temperature  rises  of  the  evolved  reaction  gases 
caused  by  exothermic  decomposition.  If  the  above  hypothesis  is 
true  then  ignition  depends  upon: 

•  temperature,  composition  and  flow  rates  of  the 
gases  evolved  by  the  heated  propellant,  and 

•  physical  and  thermal  environment  into  which  the 
gaseous  products  are  discharged. 

Fortunately,  as  will  be  seen  later  in  this  section,  one  need  not 
analyze  each  of  the  above  phenomena  to  approximate  ignition 
times.  Nevertheless  one  should  be  aware  of  the  phenomena  in 
designing  experiments  with  which  to  establish  ignition  criteria 
for  the  propellant/ conditions  of  interest. 

At  present,  ignition  criteria  are  based  upon  propellant 
temperature,  or  combinations  of  a  constant  incident  flux  and 
duration  [9].  The  former  criterion  is  of  greatest  value  in  that 
it  applies  to  time-dependent  heat  fluxes.  Ignition  is  predicted 
when  the  temperature  of  the  propellants  surface  or  at  some  depth 
exceeds  a  critical  temperature. 

In  the  remainder  of  this  section  we  shall  examine  the  effect 
of  disregarding  gas  evolution  upon  the  ignition  time.  Before 
doing  so  it  should  be  noted  that  a  one-to-one  correspondence 
does  not  exist  between  foam  temperature  and  rates  of  gas 
evolution. 
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Table  1  presents  predicted  foam  temperatures  following  ex¬ 
posure  of  HMX  to  specified  incident  heat  fluxes.  Internal  heat¬ 
ing  is  also  provided  for  even  though  it  is  not  specified  in  the 
table.  The  initial  temperature  of  the  propellant  is  294  K. 

Foam  or  propellant  surface  temperatures  are  specified  at  times 

* 

at  which  the  gas  evolution  rates  equal  specific  fractions  6  of 
the  steady  burning  rate. of  0.03  cm/sec. 


TABLE  1 .  PREDICTED  FOAM  TEMPERATURES  AND  REGRESSION  RATES  OF  HEATED  HMX 


Incident  Heat 
Flux*, 
cal/cm-sec 

Temperature  Tf  corresponding  to  specified  fraction  6 
of  the  steady  burning  rate,  K** 

6  =  0  5  =  0.1  6  =  0.5  6  =  0.9  6  -  1.0 

0.5 

555.0 

(211.44) 

565.6 

(237.19) 

580.9 

(239.33) 

587.7 

(239.58) 

588.9 

(239.60) 

1.0 

555.0 

(52.86) 

570.30 

(£1-99) 

585.9 

(64.45) 

589.1 

(64.71) 

592.3 

(64.74) 

2.0 

555.0 

(13.25) 

579.0 

(16.50) 

588.9 

(18.25) 

594.9 

(18.50) 

596.1 

(18.53) 

10.0 

555.0 

(0.529) 

603.5 

(0.668) 

607.8 

(0.875) 

611.6 

(0.959) 

612.4 

(0.974) 

30.0 

555.0 

(0.0587 

629.1 

(0.0697) 

630.4 

(0.0902) 

631.8 

(0.1031) 

632.2 

(0.1057) 

100.0 

555.0 

(0.00529) 

665.6 

(0.00577) 

664.5 

(0.00681) 

664.5  664.6 

(0.00760)  (0.00770) 

*  Incident  heat  fluxes  are  from  some  external  heat  source. 
**  Times  of  occurrence  given  in  parentheses  in  seconds. 


From  Table  1  it  may  be  observed  that  an  ignition  criterion 
based  solely  upon  a  given  temperature  implies  that  ignition 
occurs  with  differing  rates  of  gas  release.  The  latter  reempha¬ 
sises  the  question  raised  earlier  regarding  the  significance  of 
the  differing  rates  of  gas  evolution  upon  the  ignition  time. 

In  order  to  perform  the  above  assessment  we  shall  assume 
that  the  critical  ignition  temperature  lies  between  the  melt 
temperature  of  555  K  and  the  lowest  temperature  cited  in  Table  1 
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for  6  =  1.  A  6  value  of  1  Is  chosen  since  it  corresponds  to  a 
regression  rate  equal  to  the  steady  burning  rate.  Based  upon 
the  above  assumptions,  the  critical  temperature  lies  between 
555.0  and  588.9  K. 

Applying  the  above  temperature  range  to  Table  1  indicates 
that  the  ignition  time  is  between  211.4  and  239.6  sec  when  the 
incident  flux  is  0.5  Val/cm^-sec;  and  between  0.00529  and  a 
linearly  interpolated  value  of  0.00544  sec  when  the  incident 
flux  is  100  cal/cm^-sec.  Notice  that  the  above  times  differ 
by  only  13.3  and  2.8  percent  for  the  two  fluxes.  These  results 
suggest  that  ignition  temperature  provides  a  reasonably  accurate 
means  for  predicting  ignition  at  least  for  HMX.  Rates  of  gas 
evolution  are  of  secondary  importance. 

For  the  present  we  shall  assume  that  the  ignition  tempera¬ 
ture  of  HMX  equals  the  mean  temperature  (572  K)  associated  with 
the  range  of  temperatures  cited  earlier.  Ignition  times  t^  re¬ 
quired  to  achieve  a  foam  temperature  of  572  K  may  be  determined 
from  Table  1  as  a  function  of  q.  They  are  given  approximately  by 

tj.  -  59/q2  (11) 

where  t^  is  in  seconds  when  q  is  in  cal/cm  -sec.  Experimental 
results  presented  in  reference  [9]  indicate  that  the  exponent  2 
of  equation  (11)  can  vary  from  about  1.6  to  2.0  depending  upon 
the  propellant.  Most  propellants  involve  coefficients  only  a 
fraction  of  the  coefficient  59  presented  by  equation  (11) .  The 
relatively  large  coefficient  is  attributed  to  the  high  melt 
temperature  of  HMX  and  the  fact  that  much  of  the  "ignition  time" 
is  expended  in  initiating  melting.  In  this  regard,  times  to 
initiate  melting  are  proportional  to  the  square  of  the  difference 
between  the  melt  temperature  and  the  initial  temperature  of  the 
propellant  [10]. 


12 


2.4.2  Effect  of  Sudden  Pressure  Relief  Upon  KMX  Foam  Mass 


In  the  previous  interim  report  [4],  Boggs'  experimental 
results  [11]  wer-j  used  to  deduce  the  amount  of  HMX  foam  present 
during  steady  burning  as  a  function  of  pressure.  Foam  masses 
were  then  used  to  evaluate  the  heat  transfer  coefficient  asso¬ 
ciated  with  heat  flow  from  the  foam  to  the  melt  interface. 

Foam  masses  were  deduced  by  estimating  the  thickness  and 
density  of  "frozen"  HMX  layers  presented  by  Boggs  [11] .  Freezing 
was  accomplished  by  suddenly  relieving  the  pressure  (34  and  68 
bars)  present  during  steady  burning.  Questions  regarding  the 
correspondence  between  the  mass  of  "frozen"  foam  with  the  foam 
mass  present  during  steady  burning  were  not  addressed.  In  other 
words,  how  much  does  the  foam  mass  change  following  pressure  re¬ 
lief  due  to  continued  outgassing  and  melting? 

Figure  2  presents  analytical  predictions  of  the  foam  mass 
(molten  plus  frozen)  as  a  function  of  time  following  sudden  re¬ 
lief  of  the  tiro  pressures  cited  earlier.  At  time  zero  the  foam 
mass  equals  that  of  steady  burning.  Several  hundred  milliseconds 
after  pressure  relief  the  foam  is  completely  frozen. 

It  may  be  observed  that  the  foam  mass  at  34  bars  decreases 
from  0.00254  to  0.00220  g/cm^  following  pressure  relief,  while 

2 

the  foam  mass  at  68  bars  decreases  from  0.00139  to  0.00112  g/cm  . 
Percentagewise,  the  above  mass  reductions  are  13.4  and  19.4 
percent,  respectively.  The  somewhat  greater  reduction  (percent¬ 
agewise)  of  the  foam  mass  present  at  68  bars  is  due  to  its 
higher  temperature.  The  result  is  more  rapid  mass  loss  of 
initially  smaller  amounts  of  foam. 

While  the  above  mass  changes  are  significant,  they  are 
probably  small  compared  to  errors  in  estimating  the  thickness 
and  density  of  the  "frozen"  foam.  For  this  reason,  the  foam 
masses  presented  in  the  previous  interim  report  [4]  shall  not 
be  corrected  until  better  measurements  of  the  "frozen"  foam 
mass  are  available. 


otoiro 


o 


o 


o 


6 


u 


C. 


c 


3.  CRACK  PROPAGATION  INTO  HIGH  PRESSURE/ TEMPERATURE  GAS  CAVITIES 

This  report  considers  cracks  with  a  width  and  length  of 
0.2  cm  and  10.0  cm,  respectively.  The  crack  width  was  chosen 
sufficiently  large  to  accommodate  sufficient  hot  gas  to  ignite 
'  crack  surfaces.  The  crack  length  was  chosen  sufficiently  long 
to  afford  appreciable  differences  in  the  heating  of  the  crack 
surfaces , 

Initially  the  gases  within  the  crack  are  considered  to  be 
in  a  state  of  rest  at  a  temperature  and  pressure  of  294  K  and 
1  bar,  respectively.  At  zero  time  the  crack  is  considered  to 
propagate  into  a  cavity  of  gases  whose  temperature  is  3000  K. 

Gas  pressures  within  the  cavity  were  varied.  The  downstream  end 
of  the  crack  is  considered  either  closed,  or  open  to  a  second 
cavity  containing  gases  at  the  same  pressure  and  temperature  as 
those  initially  within  the  crack.  All  gases  are  assumed  to  have 
the  same  composition. 

3 . 1  Effect  of  Erosive  Heating  Upon  Steady  Burning  Conditions 

This  section  is  concerned  with  assessing  the  effect  of  ero¬ 
sive  heating  within  cracks  upon  steady  burning  of  HMX.  Table  2 
is  presented  for  reference  purposes .  It  gives  predicted  HMX 
steady  heat  fluxes  q^,  burn  rates  r^ ,  foam  masses  and  foam 
temperatures  T^  as  a  function  of  gas  pressure  P  when  erosive 
heating  is  absent.  Figures  3,  4  and  5  indicate  the  effects  of 
erosive  heating  upon  steady  bum  rates,  foam  masses  and  foam 
temperatures  as  functions  of  the  cross -flow  gas  velocity  u  for 
three  different  pressures,  namely  1,  10  and  100  bars.  The  gas 
temperature  is  held  constant  at  3000  K,  while  the  u  value  en¬ 
compasses  gas  velocities  likely  to  be  present  within  cracks. 

Figure  3  indicates  that  erosive  heating  can  increase  the  HMX 
regression  rate  by  a  factor  as  large  as  nine.  The  result  of  more 
rapid  burning  is  decreased  foam  mass  (see  Figure  4)  and  increased 
foam  temperature  (see  Figure  5).  Of  greatest  important  to  the 
DDT  problem  are  the  reductions  of  foam  masses  shown  in  Figure  4 
with  increased  gas  velocities  u.  It  may  be  observed  that  steady 
erosive  heating  can  reduce  the  steady  foam  mass  by  as  much  as  an 
order  of  magnitude. 
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TABLE  2.  STEADY  HEAT  FLUXES  (qf),  BURN  RATES  (rf),  FOAM  MASSES  (Mf ) , 
AND  FOAM  TEMPERATURES  (Tf)  DURING  HMX  BURNING  VERSUS  PRESSURE  (P) 


P,  Bars 

qf,  cal/cm2-sec 

r^,  cm/ sec 

Mc,  g/cm2 

Tff  K 

1 

10.0 

0.030 

0.0367 

601 

2 

19.0 

0.055 

0.0247 

615 

5 

44.0 

0.120 

0.0132 

635 

10 

84.0 

0.217 

0.0075 

'653 

20 

160.0 

0.394 

0.0039 

673 

30 

233.0 

0.559 

0.0027 

686 

40 

305.0 

0.716 

0.0021 

695 

50 

376.0 

0.867 

0.0016 

702 

60 

446.0 

1.010 

0.0014 

709 

70 

515.0 

1.160 

0.0012 

714 

80 

583.0 

1.300 

0.0010 

719 

90 

652.0 

1.440 

0.0009 

723 

100 

719.0 

1.530 

0.0008 

727 

3.2  Depths 

of  Propellant  Ignited  in  "Closed  Cracks" 

"Closed"  refers  to  cracks  whose  downstream  end  is  closed. 
In  such  cracks,  gases  will  ultimately  stagnate  in  the  absence 
of  burning . 


Figure  6  presents  the  maximum  depths  of  propellant  ignited 
within  HMX  cracks  as  a  function  of  the  cavity  pressure.  Depths 
are  measured  from  the  upstream  end  of  the  crack  that  is  connected 
to  the  high-pressure  cavity.  Times  at  which  the  propellant  at 
the  maximum  depths  ignites  are  presented  at  the  top  of  Figure  6. 

As  one  would  expect,  the  maximum  depth  of  propellant  ignited 
by  the  inflow  of  hot  gases  increases  with  the  cavity  pressure. 
Ignition  times  range  from  a  few  milliseconds  with  low  cavity 
pressures,  to  a  fraction  of  a  millisecond  with  high  cavity  pres¬ 
sures.  The  maximum  amount  of  foam  generated  by  the  inflow  of 

hot  gases  occurred  with  the  lowest  cavity  pressure  considered, 

/  0 

namely  S  bars.  It  was  only  6.10“  g/cnr  which  is  roughly  an 
order  of  magnitude  less  than  those  presented  in  the  next  section 


19 


Ignition  Times 


for  "open  cracks”  exposed  to  the  same  cavity  pressure.  This 
result  suggests  that  closed  cracks  are  not  apt  to  generate  ap¬ 
preciable  foam  unless  they  are: 

(1)  Sufficiently  long,  or  wide  at  their  downstream 
end  to  prevent  pronounced  decreases  of  the  heat 
flux  while  appreciable  foam  is  being  generated, 
and/or  are  able  to 

(2)  maintain  relatively  low  gas  pressures  for  suffi¬ 
cient  time  until  near  steady- state  foam  masses 
are  achieved. 

These  possibilities  will  be  explored  following  experimental 
validation  of  the  analytical  predictions. 

3.3  Predictions  for  "Open  Cracks" 

The  downstream  end  of  "open  cracks"  is  connected  to  a  cavity 
of  gas  having  the  same  pressure  and  temperature  as  the  gas  ini¬ 
tially  within  the  crack.  With  such  cracks,  hot  gases  continue 
to  flow  through  the  crack  as  long  as  the  cavities  remain  at 
different  pressures. 

Figures  7  through  9  present  gas  velocities,  crack  widths, 
gas  pressures  and  foam  masses  versus  crack  depth  for  cavity 
pressures  of  8,  16  and  68  bars,  respectively.  The  dotted, 
dashed,  and  solid  curves  refer  to  conditions  present  at  speci¬ 
fied  times  following  sudden  exposure  of  the  crack  to  the  high- 
pressure  cavity. 

Crack  widths  are  greatest  at  the  high-pressure  cavity  end 
of  the  cracks  due  to  longer  exposures  to  elevated  gas  pressures. 
Crack  expansion  will  continue  at  least  until  stress  waves  are 
applied  to  the  cracks. 

Initially  the  foam  mass  is  greatest  at  the  high-pressure 
cavity  end  of  the  crack  as  one  would  expect.  Thereafter,  the 
greatest  foam  mass  occurs  at  increasing  crack  depths.  Ultimately 
the  foam  mass  is  greatest  at  or  near  the  downstream  end  of  the 
cracks.  The  latter  is  due  to  the  fact  that  the  heat  fluxes  at 
the  downstream  end  of  the  crack  start  out  small  and  then  gradu¬ 
ally  increase.  The  converse  is  true  at  the  upstream  end  of  the 
crack. 
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Crack  Width 


Mean  values  of  the  foam  masses  presented  in  Figures  7 
through  9  are  shown  in  Figure  10  as  a  function  of  time.  The 
mean  value  equals  the  integral  of  the  foam  masses  with  respect 
to  crack  depth  divided  by  the  total  length  of  the  crack. 

Figure  10  indicates  that  the  mean  foam  mass  asymototically  ap¬ 
proaches  values  of  about  0.0040,  0.0020,  and  0.0005  g/cm^ 'for. 
cavity  pressures  of  8,  16,  and  68  bars,  respectively.  Notice 
that  the  above  foam  masses  are  roughly  inversely  proportional 
to  the  cavity  pressure.  The  greater  the  optimal  amount  of  foam 
is  the  longer  it  takes  to  generate. 

At  present  the  masses  cit^d  above  do  not  appear  adequate  to 
cause  DDT  unless  extremely  intense  stress  waves  are  generated, 
say  by  the  presence  of  multiple  burning  cracks.  Two  possibili¬ 
ties  exist  for  foam  enhancement.  The  first  is  indicated  above, 
namely  low  gas  pressures  within  the  cavity  which  lessen  the 
rate  at  which  the  propellant  is  heated.  The  second  possibility 
is  the  accumulation  of  melt  within  cracks  brought  about  by  melt 
flow  caused  by  high  velocity  gases  moving  laterally  over  the 
foam.  The  latter  possibility  remains  to  be  explored.  It  applies 
to  "closed"  as  well  as  "open  cracks". 


0.0005 


0.0002 


Time  Following  Crack  Exposure  to  Cavity  Gases,  ms 


Cavity  Gases  Initially  at  16  bars 


Time  Following  Crack  Exposure  to  Cavity  Gases,  ms 


Cavity  Gases  Initially  at  8  bars 


Time  Following  Crack  Exposure  to  Cavity  Gases,  ms 
Figure  10.  Mean  HMX  Foam  Mass  Within  "Open  Cracks"  versus  Time 
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4.  EXPERIMENTAL  DDT  DESIGN 


4.1  Design  of  DDT  Apparatus 

A  preliminary  design  of  a  fixture  for  studying  DDT  was 
presented  in  the  previous  interim  report  [4] .  Since  then  the 
fixture  has  been  revised  as  shown  in  Figure  11.  The  principal 
difference  between  the  two  designs  is  in  the  means  for  heating 
and  igniting  the  propellant.  In  the  previous  design  hot  nitrogen 
gas  flow  was  considered  as  a  heat  source  in  order  to  accentuate 
melt  formation.  Iiot  gas  flow  will  dissipate  gases  evolved  by  the 
heated  propellant  and  increase  melt  formation  by  delaying  igni¬ 
tion.  Ignition  may  be  initiated  by  turning  the  nitrogen  gas 
flow  off  after  ignition  conditions  are  achieved. 

Since  then  it  was  found  that  the  foam  mass  is  extremely 

dynamic  following  ignition  as  shown  in  Figure  12.  The  dashed 

portion  of  the  curve  shows  how  the  foam  mass  increases  while  it 

2 

is  being  heated  by  an  incident  heat  flux  of  1  cal/cm  -sec.  Ig¬ 
nition  is  assumed  to  occur  when  the  foam  mass  achieves  a  value 
almost  twice  that  of  steady  burning. 

Figure  12  indicates  that  the  foam  mass  will  rapidly  peak 
and  then  fall  below  its  steady-state  value.  Thereafter,  the 
foam  mass  will  oscillate  with  decreasing  amplitudes  about  its 
steady  value.  It  may  be  observed  that  a  time  period  of  only 
0.26  sec  elapses  following  ignition  before  the  foam  mass  de¬ 
creases  to  its  steady  value .  Initiating  closure  of  the  gas 
space  within  such  a  short  time  period  raises  instrumentation 
difficulties.  The  first  is  in  the  time  expended  in  detecting 
ignition  and  in  mechanically  releasing  the  driver.  These  events 
would  require  about  0.1  sec.  Secondly,  a  delay  of  about  0.2  sec 
would  be  required  before  the  driver  strikes  the  piston.  While 
one  can  substantially  reduce  the  free-fall  time  of  the  driver, 
say  by  the  use  of  compressed  gas,  it  is  not  warranted  at  the 
present  time.  In  this  regard,  initial  tests  should  be  kept  as 
simple  as  possible  until  the  design  functions  as  anticipated 
and  DDT  is  achieved. 
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Figure  11.  Schematic  of  Revised  DDT  Apparatus 


Time  Following  Start  of  Heating,  sec 


For  the  above  reasons  it  is  planned  to  burn  the  propellant 
steadily  before  releasing  the  driver.  While  steady  burning  does 
not  provide  for  optimum  foam  masses,  it  has  two  benefits.  The 
first  is  consistent  burning  conditions  at  the  start  of  closure 
of  the  gas  space.  The  second  is  in  being  able  to  use  a  much 
simpler  heat  source,  namely  the  electrical  heater  illustrated 
in  Figure  11. 

Immediately  following  ignition  the  heater  must  be  withdrawn 
and  the  burning  allowed  to  proceed  until  steady  burning  is 
achieved.  Thereafter  the  driver  is  released  forcing  the  piston 
dovTrward.  After  passing  the  gas  vent  the  piston  will  commence 
to  increase  the  pressure  within  the  gas  space.  Closure  rates 
may  be  varied  by  altering  the  initial  heights  h^Q  and  h  .  By 
such  means  it  is  possible  to  approximate  closure  rates  of  cracks 
caused  by  stress  waves  with  amplitudes  ranging  up  to  a  few 
hundred  bars . 

Pressure  rises  will  accentuate  propellant  heating  by  in¬ 
creasing  the  pressure-dependent  term  of  equation  (10).  The  re¬ 
sult  will  be  pronounced  increases  in  the  rates  of  burning  and 
pressure  rises  that  will  continue  to  accelerate  until  the  foam 
is  essentially  consumed.  While  gas  pressures  are  rising,  stress 
waves  will  radiate  from  the  gas  space.  Following  reflection 
from  the  steel  base  beneath  the  propellant,  the  stress  waves  will 
return  and  commence  to  compress  the  gas  space,  Transit  times 
of  the  radiating  and  reflected  stress  waves  are  of  the  order  of 
10  ysec.  They  may  be  varied  by  adjusting  the  height  hpQ  of  the 
propellant  cylinder. 

4.2  Analytical  Results 

Five  factors  affect  pressure  buildup  within  the  gas  space. 
They  are : 

(1)  Initial  altitude  (h ,  )  of  the  driver  above  the 

piston,  Q0 

(2)  Initial  height  (h  )  of  the  gas  space, 

(3)  Initial  height  (h  )  of  the  propellant  cylinder, 


(4)  Initial  foam  mass  (M^) , 

(5)  Initial  foam  temperature  (T^),  and 

(6)  Driver  mass  (W) 

The  remainder  of  this  section  is  concerned  with  the  conse¬ 
quence  of  varying  items  1  through  4  individually.  Nominal  values 
for  the  above  parameters  are  presented  below: 

hdo  ■  1  ft 

h  =  2.5  cm 

go 

h _ =  5.0  cm 

po 

Mf  =  0.0367  g/cm2 

Tf  =  601  K 

W  =  100  lb 

The  above  combination  of  parametric  values  was  chosen  in 
that  it  results  in  near-optimal  pressure  rises  in  the  gas  space. 
Values  assigned  for  the  mass  and  temperature  T^q  are  asso¬ 
ciated  with  the  steady  burning  of  HMX  at  1  bar  of  pressure  (see 
Table  2). 

Figure  13  presents  predicted  gas  pressures  P  and  the  gas 
space  heights  h^  as  a  function  of  time  until  the  pressure  peaks. 
Figure  14  presents  extremely  dynamic  portions  of  the  Figure  13 
curves  on  either  side  of  the  pressure  peak. 

Figure  13  indicates  the  pressure  will  peak  at  about  30  kbars 
in  10  msec  following  the  start  of  closure  of  the  gas  space. 

Figure  14  suggests  that  much  of  the  above  pressure  rise  will 
occur  within  only  a  fraction  of  a  ysec.  Thereafter  the  pressure 
decays  in  an  exponential  fashion  as  a  consequence  of  expansion 
of  the  gas  space  by  the  pronounced  gas  pressures. 


Results  shown  in  Figures  13  and  14  are  predicated  on  no  heat 
losses  from  the  foam  to  the  piston,  and  the  presence  of  planar 
stress  waves  whose  amplitude  is  not  altered  in  moving  through 
the  propellant.  There  is  a  distinct  possibility  of  foam-piston 
contact  in  that  the  surface  of  the  burning  propellant  will  become 


concave.  Moreover  the  foam  layer  is  inherently  irregular.  It 
is  estimated  that  the  foam  may  be  in  contact  with  the  pi3ton  for 
as  long  as  1  ms  and  thereby  lose  a  significant  fraction  of  its 
sensible  heat.  As  a  consequence  the  bottom  face  of  the  piston 
should  be  covered  with  a  thin  layer  of  high  temperature  insula¬ 
tion.  Planar  stress  waves  may  be  approximated  by  using  propel¬ 
lant  cylinders  with  relatively  large  cross-sectional  areas  for 
the  given  cylinder  length. 

Figures  15  through  17  illustrate  the  consequence  of  varying 
the  initial  values  for  h^,  h  ,  and  h^  upon  the  peak  gas  pressure 
and  associated  void  height.  Pressure-time  curves  are  not  pre¬ 
sented  in  that  their  shape  is  very  similar  to  that  shown  in 
Figure  14.  In  the  above  figures  plots  of  the  maximum  pressure 
achieved  and  corresponding  height  of  the  gas  space  are  given  as 
functions  of  the  initial  value  of  the  specified  parameter  under 
consideration.  Remaining  parameters  are  equal  to  the  nominal 
values  cited  earlier. 

Figure  15  indicates  that  the  driver  should  be  at  least 
0.5  ft  above  the  piston  in  order  to  achieve  near  optimal  gas 
pressures.  Peak  pressures  do  not  differ  appreciably  for  driver 
altitudes  between  about  0.5  and  2  ft. 

Figure  16  indicates  that  optimal  value  of  the  peak  pressure 
will  be  produced  when  initial  gas  space  height  h^o  is  in  the 
vicinity  of  2.5  cm.  Above  and  below  2.5  cm,  the  pressure  de¬ 
creases  monotonically .  In  this  regard  smaller  heights  h^o 
lessen  the  time  period  over  which  the  gas  space  is  compressed. 
The  result  is  lower  foam  temperatures  and  hence  less  dynamic 
burning.  On  the  other  hand,  heights  hgQ  greater  than  2.5  cm  in¬ 
crease  the  gas  volume  at  the  time  at  which  the  peak  pressures 
occur,  and  thereby  lessen  the  peak  pressure. 

Figure  17  illustrates  the  consequence  of  varying  the  ini¬ 
tial  height  hpQ  of  the  propellant  cylinder.  It  may  be  observed 
that  decreases  in  the  height  of  the  propellant  cylinder  accentu¬ 
ate  the  gas  pressures.  Increased  gas  pressures  are  caused  by 
more  rapid  reflection  and  return  of  stress  waves  radiated  from 
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the  gas  space.  The  result  is  higher  amplitude  stress  wave's 
being  applied  to  the  gas  space  at  least  during  periods  of  rising 
pressure. 

Figure  18  shows  the  effect  of  varying  the  initial  foam  mass 
M^.  It  may  be  observed  that  the  peak  pressure  is  highly  depen¬ 
dent  upon  the  foam  mass.  That  is  why  we  are  so  concerned  with 
the  generation,  and  possible  transport  and  accumulation  of  malt 
within  cracks . 

Finally  it  is  important  to  know  how  well  the  initial  tempera 
ture  of  the  propellant  needs  to  be  controlled  to  achieve  consis¬ 
tent  test  results.  Of  concern  is  the  sensitivity  of  the  steady 
bum  rate,  and  mass  and  temperature  of  the  foam  mass  to  the 
initial  temperature  of  the  propellant.  Figure  19  illustrates 
the  above  dependence.  Notice  that  the  foam  mass  is  optimal 
with  initial  propellant  temperatures  between  about  270  and  320  K. 
High  initial  temperatures  accentuate  the  burn  rate  as  well  as 
the  foam  temperature.  For  purposes  of  test  consistency,  the 
initial  temperature  of  the  propellant  should  be  maintained  within 
about  +20  K  of  the  normal  room  temperature  of  294  K.  Only  during 
very  cold  days  will  temperature  control  be  needed. 
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5.  SUMMARY ,  CONCLUSIONS,  AND  FUTURE  PLANS 


5.1  Crack  Studies 

These  studies  were  concerned  with  assessing  the  consequence 
of  crack  propagation  into  cavities  containing  gases  at  higher 
pressures  and  temperatures  than  initially  exist  within  the  crack?. 
The  result  of  this  and  previous  studies  [2-4]  reveal  three  im¬ 
portant  factors  in  causing  DDT.  They  are  the  presence  of: 

(1)  appreciable  amounts  of  foam  or  melt,-" 

(2)  high  temperature  foam,  and 

/ 

(3)  high  amplitude  stress  waves,..  ' 

Item  1  is  important  in  that  much  of  the  foam  is  consumed 
during  accelerated  burning.  The  latter  includes  the  foam  present 
at  the  start  of  the  pressure  transient  plus  additional  foam 
generated  by  melting  during  the  pressure  transient.  Usually  the  “ 
added  foam  is  a  small  fraction  of  that  present  at  the  start  of 
the  pressure  transient. 

Item  2  is  important  in  that  high  foam  temperatures  speed 
foam  consumption  and  thereby  decrease  the  time  over  which  cracks 
are  expanded  by  elevated  gas  pressures.  The  importance  of  item  3 
lies  in  the  fact  that  high  amplitude  stress  waves  not  only  trigger 
accelerated  burning  but  also  retard  crack  expansion.  In  this 
regard,  gas  pressures  always  peak  while  cracks  are  expanding 
rather  than  contracting. 

Analytical  studies  reveal  that  single  burning  cracks  are 
more  likely  to  cause  DDT  when  they  are  in  close  proximity  to 
the  motor  case.  Close  proximity  leads  to  more  rapid  reflection 
and  return  of  stress  waves  radiating  from  the  crack.  The  result 
is  more  intense  stress  waves  being  applied  to  the  crack  at  least 
during  periods  of  rising  pressure.  In  turn,  high  amplitude  stress 
waves  accentuate  gas  pressures  within  the  cracks  by  speeding  the 
burning  and  restricting  crack  expansion. 


The  above  DDT  possibility  involves  only  a  single  crack.  It* 
is  in  addition  to  the  DDT  possibility  advanced  earlier  [4]  for 
situations  in  which  multiple  burning  cracks  are  present.  In  the 
latter  situations,  stress  waves  can  be  progressively  amplified 
in  moving  from  one  crack  to  the  next  crack  and  thereby  progres¬ 
sively  enhance  gas  pressures  produced  within  the  cracks. 

Several  propellant  properties  affect  the  pressure  transients 
produced  in  burning  cracks.  The  first  is  the  heat  transfer  co¬ 
efficient  associated  with  the  transfer  of  heat  from  the  foam  to 
the  solid  propellant  interface.  Greater  heat  transfer  coeffi¬ 
cients  increase  foam  buildup.  Of  equal  importance  are  the 
propellants  activation  energy  and  internal  heat  constant.  In 
this  regard,  large  activation  energies  or  internal  heats  act  to 
accelerate  the  burning. 

At  present  we  are  primarily  concerned  with  phenomena  affect¬ 
ing  foam  generation  within  cracks.  For  this  purpose  we  have  ex¬ 
amined  the  generation  of  foam  within  cracks  that  propagate  into 
cavities  containing  gases  at  much  higher  temperatures  and  pres¬ 
sures  than  initially  exist  within  the  cracks.  The  downstream  end 
of  the  cracks  was  considered  connected  to  a  low  pressure  cavity 
so  that  gases  continued  to  flow  through  the  crack.  Not  surpris¬ 
ingly,  foam  masses  decreased  with  increased  pressures  within  the 
cavity  at  the  upstream  end  of  the  crack.  Thus,  cavity  pressures 
of  8,  16  and  68  bars  yielded  average  foam  masses  of  approximately 
0.004,  0.002,  and  0.0005  g/cm^,  respectively.  Notice  that  the 
above  foam  masses  are  roughly  inversely  proportional  to  the 
cavity  pressure. 

In  the  previous  interim  report  [4]  it  was  shown  that  pres¬ 
sures  of  the  order  of  tens  of  kbars  may  be  generated  within 
multiple  burning  cracks  contained  0.005,  0.010,  and  0.015  g/cm^ 
of  foam.  Less  cracks  are  needed  to  develop  such  pressures  the 
greater  the  amount  of  foam  present.  Most  of  the  pressure  rise 
occurred  within  1  ysec  or  less. 


5.2  Peripheral  Studies 


Additional  studies  involved  predicting  the  ignition  of  HMX, 
errors  in  estimating  foam  masses  by  pressure  relief  quenching 
of  the  burning,  and  the  consequence  of  erosive  heating  upon  foam 
masses. 

Cursory  examination  of  HMX  ignition  suggests  that  foam 
temperature  presents  a  satisfactory  criterion  for  predicting  ig¬ 
nition.  Analysis  indicates  that  the  ignition  time  of  HMX  is 
inversely  proportional  to  the  square  of  the  incident  heat  flux. 
Results  remain  to  be  confirmed  experimentally. 

Pressure  relief  has  a  small  but  significant  effect  upon  HMX 
foam  mass  while  it  undergoes  solidification.  Analysis  indicated 
that  the  foam  mass  is  reduced  by  13  and  19  percent  during  pres¬ 
sure  relief  extinguishment  of  steady  burning  at  34  and  68  bars, 
respectively. 

Erosive  heating  is  due  to  the  interaction  of  gases  leaving 
a  propellant  with  gases  flowing  laterally  through  a  crack.  It 
was  found  that  erosive  heating  very  appreciably  reduces  the 
generation  of  foam.  HMX  foam  masses  can  be  decreased  by  as  much 
as  one-tenth  of  their  normal  burning  (no  lateral  gas  flow)  values 
depending  upon  the  rates  of  mass  transport  by  the  two  gas  flows. 

5.3  Analysis  of  Experimental  DDT  Design 

An  experimental  fixture  for  studying  DDT  was  described  (see 
Figure  11).  Assessments  were  made  of  the  influence  of  various 
fixture  parameters  upon  pressure  transients  produced  by  com¬ 
pressing  the  gas  space  over  burning  HMX.  It  was  found  that 
pressure  rises  of  the  order  of  a  few  tens  of  kbars  may  be  gener¬ 
ated  within  fractions  of  1  ysec.  DDT  appears  likely  with  such 
rapid  and  pronounced  pressure  buildups  [15]. 

5.4  Future  Plans /Needs 

During  the  forthcoming  year  it  is  planned  to: 

(1)  conduct  the  DDT  experiments  described  in  Section  4 
using  a  composite  propellant  consisting  of  65  per¬ 
cent  HMX  and  35  percent  polyethyleneglycol 
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(2)  analyze  possibilities  by  which  melt  may  accumu¬ 
late  within  cracks  due  to  possible  displacement 
of  the  melt  by  gases  flowing  along  the  length 
of  the  crack 

(3)  measure  foam  masses  produced  by  steady  burning 
of  HMX  and/or  the  composite  propellant  described 
by  item  1  at  1  atmosphere  of  pressure. 

The  initial  DDT  experiments  (item  1)  are  primarily  concerned 
with  establishing  whether  or  not  DDT  occurs.  Expensive  pressure 
transducers  shall  be  reserved  for  follow-on  experiments  if  DDT 
occurs.  Item  2  involves  development  of  a  subroutine  for  the 
computerized  model  PROS  described  in  sections  1  and  2.  Item  3 
is  important  in  establishing  whether  ot  not  extrapolated  HMX  foam 
masses  used  in  this  report  are  reasonable.  This  comparison  is  also 
needed  in  that  the  heat  transfer  coefficient  associated  with  heat 
flow  from  the  foam  to  the  melt  interface  was  deduced  from  the 
pressure  dependent  foam  mass  present  during  steady  burning. 

Follow-on  plans  will,  of  course,  depend  upon  information 
gathered  from  the  above  activities.  At  a  minimum  they  should 
include : 

•  upgrading  DDT  experimental  fixture  (see  Figure  11) 
to  provide  for  variations  of  the  initial  foam  mass. 

This  may  be  achieved  by  controlling  the  size  of 
the  gas  vent  and  hence  the  pressure  in  which  the 
propellant  bums  steadily  prior  to  releasing  the 
driver. 

•  conducting  DDT  experiments  with  various  composite 
propellants  containing  HMX 

•  using  analysis  and  experimental  results  to  iden¬ 
tify  compositions  and  critical  property  values 
of  composite  propellants  that  should  be  avoided 
to  minimize  the  likelihood  of  DDT. 

Follow-on  experiments  should  be  primarily  concerned  with  estab¬ 
lishing  threshold  "crack",  burning  and  motor  conditions  needed 
to  initiate  various  propellants.  Analysis  is  needed  to  estab¬ 
lish  the  likelihood  of  achieving  the  threshold  conditions. 
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APPENDIX:  NOMENCLATURE  AND  HMX  DATA 


This  appendix  describes  variables  used  in  this  report  as 
well  as  the  properties  used  for  HMX  and  its  reaction  gases.  The 
latter  is  presented  in  Table  3. 


Nomenclature 

a  rate  of  steady  burning  at  1  bar  of  pressure,  cm/sec 

C  specific  heat  of  gases  evolved  by  propellant  at  constant 

s  pressure,  cal/g-  K 

Cm  specific  heat  of  molten  propellant,  cal/g-  K 

Cp  specific  heat  of  solid  propellant,  cal/g-  K 

C  crack  width,  cm 

w 

c,  ,c.?  constants  used  to  describe  heat  transfer  coefficient  h 
where  h  -  c^(Z  exp  -  E/Tf)c2 

E  activation  energy  of  propellant  divided  by  gas  constant, 

K 

h  heat  transfer  coefficient  associated  with  heat  flux  from 

foam  to  melt  interface,  cal/cm^-sec-  K 

2 

hc  convective  heat  transfer  coefficient,  cal/cm  -sec-  K 
h^Q  initial  altitude  of  driver  before  piston  (see  Figure  11) , 


k 


height  of  gas  space,  cm 

initial  value  of  h  (see  Figure  11) ,  cm 

D 

height  of  propellant  cylinder,  cm 

initial  value  of  h^  (see  Figure  11) ,  cm 

mechanical  impedence  of  solid  propellant,  equals  I  + 
0.0002  P,  bars-sec/cm 

see  I 

constant  used  to  predict  erosive  burning  rate  in  equation 
(4),  cm3-  K/cal 


thermal  conductivity  of  solid  propellant,  cal/cm-sec-  K 
indepth  crack  distance,  cm 

2 

mass  of  unit  area  of  foam  layer,  g/cm 

molecular  weight  of  reaction  gases ,  g/mole 

exponent  of  pressure  P  used  to  describe  steady  burning 
rate  r£ 

gas  pressure,  bars  or  kbars 
maximum  value  of  P  achieved,  kbars 
Prandtl  number,  dimensionless 

sensible  heat  per  unit  mass  of  foam  during  steady  burning, 
cal/g 

latent  heat  of  fusion  of  propellant,  cal/g 

heat  of  reaction  of  propellant  gases  cal/g 

heat  generated  within  foam  per  unit  mass  of  evolved  gas, 
cal/g 

heat  flux  incident  upon  propellant  from  external  heat 
sources,  cal/cm2-sec 

rate  of  heating  of  unit  area  of  foam,  eal/cm2-sec 

2 

steady  value  of  q^,  cal/cm  -sec 
rate  of  propellant  burning,  cm/ sec 
steady  value  of  r^  at  given  pressure,  cm/ sec 
gas  constant,  cm/  K 

temperature  of  foam  or  propellant  surface,  K 

foam  temperature  during  steady  burning  at  given  pressure, 
K 

temperature  of  combustion  gases,  K 
melt  temperature  of  propellant,  K 
initial  temperature  of  propellant,  K 
time,  sec 


ignition  time,  sec 

u  cross  flow  gas  velocity,  cm/ms  or  cm/ sec 
W  mass  of  driver  (see  Figure  11),  lb 
Z  frequency  factor,  1/sec 

6  constant  used  to  describe  erosive  heating  (see  equation 

(10)),  dimensionless 

3 

d  density  of  solid  propellant,  g/cm 

y  ratio  of  specific  heats  of  combustion  gases,  dimensionless - 

TABLE  3.  PROPERTIES  OF  HMX  PROPELLANT  AND  GASES 


Parameter 

Value 

Source 

a,  n 

a  =  0.030  cm/sec 
n  «=  0.86 

Reference  [12] 
Reference  [12] 

cl’c2 

ci  =  7.6  •  10^  cal/cm^-sec-  R 

C2  =  0.338  (dimensionless) 

Calculated  [4] 
Calculated  [4] 

C 

P 

0.4  cal/g-  K 

Reference  [13] 

0.5  cal/g-  K 

Reference  [13] 

E 

27,000  K 

Reference  [14] 

I 

o 

0.450  bars-sec/cm 

Assumed 

K 

0.0013  cal/cra-sec-  K 

Assumed 

M 

w 

24.1  g/mole 

Assumed 

P 

r 

0.7  (dimensionless) 

Calculated 

Sn 

50  cal/g 

Reference  [13] 

Qr 

1300  cal/g 

Assumed 

Qs 

150  cal/g 

Assumed 

R 

3517  cm/  K 

Assumed 

ni 

555  K 

Reference  [13] 

z 

0.5  •  1020/sec 

Reference  [14] 

6 

15  (dimensionless) 

Assumed 

P 

1.9  g/cm2 

Reference  [12] 

Y 

1.2  (dimensionless) 

Assumed 
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